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ABSTRACT 

The process of turbulent radial mixing in protoplanetary disks has strong relevance to the analysis of the spatial 
distribution of crystalline dust species in disks around young stars and to studies of the composition of meteorites and 
comets in our own solar system. A debate has gone on in the recent literature on the ratio of the effective viscosity 
coefficient i/ (responsible for accretion) to the turbulent diffusion coefficient D (responsible for mixing). Numerical 
magneto- hydrodynamic simulati ons ha ve yielded values between f/D ~ 10 (Carballido, Stone & Pringle [20051 1 and 
u/D ~ 0.85 (Johansen & Klahr [20051 ) . Here we present two analytic arguments for the ratio v/D = 1/3 which are 
based on elegant, though strongly simplified assumptions. We argue that whichever of these numbers comes closest 
to reality may be determined observationally by using spatially resolved mid-infrared measurements of protoplanetary 
disks around Herbig stars. If meridional flows are present in the disk, then we expect less abundance of crystalline dust 
in the surface layers, a prediction which can likewise be observationally tested with mid-infrared interferometers. 
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1. Introduction 

Turbulent mixing plays an important role in the physics of 
protoplanetary accretion disks. The same turbulence that is 
responsible for the anomalous viscosity of the disk (and thus 
for the accretion process) is also responsible for the radial 
and vertical mixing of material in the disk. This mixing is 
likely to have strong influence on thermochemical processes 
and grain growth in disks. For instance, models of the gas- 
phase chemistry in disks turn out to be strongly dependent 
on the level of turbulent mixi ng (G ammic 199G; Semenov, 
Wiebe & Henning(2003; Ilgner [20061 ) . This is because some 
chemical reactions may be very slow in some regions but 
fast in another. Thus mixing can transport processed gas 
from the latter region to the former region, affecting the 
molecular abundances there. 

For dust solid-state chemistry mixing plays presumably 
an even bigger role. Some evidence from meteoritics points 
to weak mixing, such as the recently discovered chemical 
complementarity of chondrules and matrix in some chon- 
drites (Klerner & Palme [20001 ). On the other hand, the 
existence of Calcium-Aluminium-rich Inclusions (CAIs) in 
many chondritic meteorites points toward some level of 
mixing in the early solar system, since CAIs are thought 
to form closer to the sun. Samples returned from comet 
Wild 2 by the STARDUST mission revealed that CAI-like 
material can be present in comets, which is interp reted as 
evidence for radial mixing (Zolensky et al. [2006( ). In in- 
frared observations of protoplanetary disks there is strong 
evidence of the existence of crystalline silicates at radii in 
the disk that are far larger than the radius at which the 
disk temperatures are high enough for thermal annealing. 
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One idea is that they got there by radial mixing from these 
hot inner region s to t he cooler outer regions (Gail [200 1[ : 
Bockelee-Morvan [20021) . 

The process of radial mixing by turbulence can be de- 
scribed by a diffusion equation governed by a diffusion co- 
efficient D (MorfiU & Volk ,1984 ). Turbulence also plays 
an important ro le in the theory of coagulation of dust 
(e.g. Volk et al. [lOMl . It can both prevent and acceler- 
ate grain growth, in complex ways. The process of accre- 
tion is, on the other hand, described by the equations of 
mass- and angular momentum conservation in a viscous 
medium, which is ruled by the e ffective viscosity coeffi- 
cien t V (Shakura & Sunvaev [1973[ : Lynden-Bell & Pringle 
1974). Both v and D may depend on distance to the star 
(i.e. the radial coordinate of the disk R) and both have the 
same dimension. In fact, because they both arise from the 
same turbulence, there are reasons to believe that they are 
nearly the same, apart from a dimensionless factor of order 
unity which is called the Schmidt number: v/ D = Sc. The 
effective viscosity v may be partly produced by Reynold 
stress (turbulent motions of the gas). Maxwell-stress (mag- 
netic field lines transporting angular momentum), as well as 
gravity waves in moderately gravitationally unstable disks. 
However, the mixing D can only be due to Reynold stress 
(turbulent motions). So D and v are clearly related, but do 
not necessarily have to be equal. 

As has been shown by Clarke & Pringle (jl988t ). the 
efficiency of outward radial mixing in steady disks de- 
pends on the ratio of v to 13, i.e. on Sc, and therefore 
the Schmidt number plays an essential role in the under- 
standing of the distribution of various chemical and dust 
species in disks, much more so than the absolute values of 
D and v. The Schmidt number of turbulent flows has been 
discussed in several studies in the past (e.g. Tennekes & 
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Lumlev [l97l McComb llQQOf) . These authors typically ar- 
gue for u/ D ~ 0.7. It is not clear, however, whether these 
findings can be directly translated to turbulence in rotat- 
ing Keplerian accretion disks. Accretion disk flows a re fun - 
damentally diffe rent f rom laboratory flow (Balbus 120031 ). 
Stone & Balbus (|1996D have shown, using 3D hydrodynam- 
ical simulations as well as analytical arguments, that hy- 
drodynamical turbulence resulting from vertical convection 
tends to transport angular momentum inward instead of 
outward. This stands in contrast to turbulence in planar 
shear flows which, as Stone & Balbus confirm in their sim- 
ulations, transports momentum in the opposite direction. 
The same kind of inward an gular momentum transport was 
reported by Riidiger et al. ( 20051 ). 

Currently the most widely accepted theory for the origin 
of effective viscosity in accretion disks is that of magneto- 
rotational turbulence. This picture is based on an interplay 
between weak magnetic fields and the differentially rotat- 
ing gas in the disk, causing a n inst ability that drives the 
turbulence (Balbus & Hawlev [l99lh . Angular momentum 
is then transported both by Reynold stresses as well as by 
Maxwell stresses. However, there are various uncertainties 
about whether the magneto-rotational instability can oper- 
ate in protoplanetary disks whi ch a re often very near to be- 
ing neutral (e.g. S emenov et al. 12004 : Ilgner & Nelson [20061 
Oishi et al. |2007[ ) . All in all it is still quite unclear what 
the nature of the turbulence and the mechanism of angular 
momentum transport is, and therefore the issue of the v/ D 
ratio also remains open. 

Recently, several publications have attempted to shed 
light on the v/D ratio in accretion disks with full 3D 
MHD simulations of such turbulence. These simulations 
have yielded values ranging from v/D ~ 10 (Carballido 
et al. |2005[ ) to v/D ~ 0.85 (Johansen & Klahr t2005i) . The 
drawback of these 3D MHD slab geometry models is that 
they have a finite numerical resolution which is in general 
rather coarse. There are reasons to believe that simulations 
with higher resolution might yield different results. Such 
initial calculations seem to indeed indicate a reduction of 
magneto-rotational instability, (Dzyurkevich 2007, private 
communication) . 

A t the same time, new 2D models of Keller & 
Gail (|2004D and Tscharnuter & Gail (2007) show the pres- 
ence of large-scale circulations within the disk. These calcu- 
lations support the results found previously in the analyti- 
cal paper of Urpin (Il984[) . and in the asymptotic study of 
Regev & Gitelman (2002). Following these results, beyond 
the certain critical radius near the disk's equatorial plane, 
the material moves in the outward direction, whereas the 
accretion flow develops in the surface layer of the disk. As 
a result of the large-scale circulation, which is driven by 
viscous angular momentum transfer, advective transport 
dominates diffusive mixing in the outer part of the disk. 
Species that are produced or undergo chemical reactions 
in the warm inner zones of the disks are advectivelly trans- 
ported into the cool outer regions. Such dynamics cannot be 
purely represented by any ID viscous or diffusional models. 
Thus, v/D ratio as direct measure of the mixing efficiency 
is not appropriate in this picture. But, despite the beauty 
of this theory, it is based on the postulated anomalous a- 
viscosity, which is an apriory micro-scale quantity. If strong 
turbulence exists on the scale which is comparable to the 
height of the disk, it may destroy such a circulation pattern. 



The question is now: which scenario of the disk evo- 
lution (viscous, diffusion or advective) is the right one in 
nature? From theoretical arguments it is still challenging 
to tell. But we propose here an observational test that may 
possibly distinguish between all these cases. This allows us 
to 'measure' effective v/D for stars as bright as Herbig Ae 
stars. Unfortunately this will not work for T Tauri stars 
and Brown Dwarfs, due to lack of spatial resolution. 

This paper is organi zed as follows. In Sect. [2] we re- 
derive the Pringle (|1981h equation on the basis of a simple 
diffusion recipe, and we thereby obtain the Schmidt number 
under the assumption that the simple diffusion recipe de- 
scribes the true motion of gas parcels in the disk. In Sect. [3] 
we derive the Schmidt number in a different way, by assum- 
ing that each component of the fluid obeys the same equa- 
tion. We obtain the same Schmidt number as in Sect. [3 
In Sect. 3] we review how the Schmidt number affects the 
abundance of crystalline silicates in the inner disk regions. 
In Sect. [5] we present radiative transfer results for the mid- 
infrared spectra of spatially resolved disks, and propose how 
the study of such spectra might provide insight into what 
is the value of the Schmidt number in real protoplanetary 
disks. 

In this paper we focus purely on the diffusion of gas or 
of small particles that are well coupled to the gas. Bigger 
dust particles (in the size range of 1 cm or bigger) will 
decouple from the turbulence and drift inward. This is an- 
other process which we do not include in this paper. Also, 
for simplicity we assume that the disk is vertically averaged 
and axially symmetric. 

2. The difFusion equation for the Keplerian disk: 
derivation from a discrete model 

In this s ection we will re-derive the well known equation of 
Pringle (|l98lh for the spreading of the accretion disk. The 
derivation is based on the idea that the evolution of the 
disk can be completely described by a special kind of tur- 
bulent motion of fluid parcels. In normal turbulence each 
eddy of the disk has a deviation 8v from the average veloc- 
ity of the fluid. These are generally non-Keplerian. But the 
conservation of angular momentum also prevents these gas 
parcels from drifting too far inward or outward, unless they 
exchange angular momentum with neighboring gas parcels, 
which leads to angular momentum transfer in the wrong di- 
rection (see above). However, it is known that weakly mag- 
netized disks are unstabl e to t he magneto-rotational insta- 
bility (Balbus & Hawley 1991). Magnetic field lines enable 
two parcels at different radii to exchange angular momen- 
tum 'over a distance', if they are thread by the same field 
lines. We suggest an extremely simplified picture for the 
motions of these parcels: Consider two parcels, one slowly 
drifting inward, the other slowly drifting outward. We as- 
sume that they both have a locally Keplerian angular ve- 
locity at all times. This means that the inward moving par- 
cel loses angular momentum while the other gains angu- 
lar momentum. We assume that neither parcel has angu- 
lar momentum exchange with its surroundings. So angular 
momentum can only be conserved if the loss of angular mo- 
mentum of the inward moving parcel is compensated by the 
gain of angular momentum by the outward moving parcel. 
We invoke an unspecified long-range force that enables this 
exchange of angular momentum between the two parcels. 
Physically, the only long-range force that might facilitate 
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this is the Lorenz force, by magnetic fields, so we assume 
that magnetic fields indeed provide this torque. The accre- 
tion and spreading of the disk is, in this picture, nothing 
other than a material diffusion: Packages of gas being trans- 
ferred randomly, under the condition of local angular mo- 
mentum conservation. The radial mixing of a passive tracer 
in the disk is therefore an integral part of this. If the accre- 
tion/spreading of the disk is then viewed in a framework 
of 'viscous disks', then the v/D ratio follows strictly from 
this derivation. 

This pictur e wa s also presented in Tutukov & 
Pavlyuchenkov (j2004f ). where different types of astrophys- 
ical disks are described in frames of non-stationary diffu- 
sional models. However, they used only numerical models 
and did not provide the equation which governs the evolu- 
tion of such disks. 

In this section we will derive these equations and cast 
them in the standard form of viscous disk equations, which 
yields the ratio v/D. We are fully aware that the simplifica- 
tions of our assumptions are quite drastic and that realistic 
MHD simulations are presumably better and yield some- 
what different results. Nevertheless, we justify our approach 
in light of the fundamental nature of these equations. 

2.1. The mass fluxes from the cell 

Let the mass of the disk be small enough so its self- 
gravitation can be neglected. The disk is assumed to be 
geometrically thin, Keplerian, and axially symmetric. We 
describe the structure of the disk in terms of the surface 
density S(-R), where R is the distance from the axis of sym- 
metry. 

We introduce a grid {Ri}, which divides the domain 
into annular elements which form the grid cells. We assume 
that the cell size hi={Ri-\-i-Ri) is equal to the characteristic 
radial mean free path of a turbulent element. The surface 
density is assumed to be constant within each cell. Along 
with Ej, the average orbital Keplerian velocity is also de- 
termined in each grid cell, Vi^ as well as the average radial 
component of the turbulent velocity, 14-,turb(*)- 

Let us consider a current cell and denote the nearest 
cells A and B, see Fig. [TJ We suppose that in a time At 
the mass M in the current cell moves to the neighboring 
cells due to turbulent motion. The matter leaving this cell 
is transferred to the two adjacent cells, and this redistri- 
bution of the mass must conserve mass and angular mo- 
mentum according to the philosophy described above. The 
corresponding equations for the system of three cells can 
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Fig. 1. The scheme of the mass transfer in the model of 
diffusion disk. 



M3(R) 



R R+h 



M,(R+h) 

Fig. 2. The scheme of the mass transfer between two cells 
in the model of diffusion disk 



be written 



M = Ma + Mb 

MRV = MaRaVa + MbRiVb, 



(1) 

(2) 



where Va and Vb are the Keplerian velocities for radii Ra 
and Rb, respectively. 

Here we will derive the differential equation which de- 
scribes the evolution of such a disk in the limit h ^ 0, 
assuming also for simplicity that /ij=/i=const. Using the 
Taylor expansion of Eq. ^ we obtain (see Appendix Al): 



M 



1 + 

4i? 



1 ^ 

Aff, = — 1 

2 \ 4i? 



(3) 
(4) 



Note that this looks as if mass is always transported out- 
ward, since Alb < Ma- However, as we shall see below, the 
true radial flux of matter may have either sign. 

2.2. The mass flux across the cell boundary 

Denote M^^^ as a mass transferred from the cell with R 
into the cell R + h. Also, let M^^^''^ be a mass transferred 
from the cell R + h into the cell R, see Fig. [2] Denote also 
Fab as a total flux across the boundary between the cells: 



FafcAt = m(^) - Aff +''^ (5) 

Using the Taylor expansion of Eq. 
we flnd: 

FabAt^-^hR'^'-^iMR-'/'). 

Note, that if the mass flux M'--^^ does not depend on radius, 
then Fab > 0, i.e. the mass is transferred outwards. 

One of the basic assumptions of our diffusion model is 
the relation for the flux from the cell: 



(see Appendix A2) 
(6) 



M = 2ttR • E • V; turb • At. 



(7) 



If we substitute this relation into Eq. ([6]) and introduce 
the diffusion coefficient (sec Appendix A3 for the choice of 
numerical coefficient): 



D 



1 



-h-Vr 



turb ; 



(8) 
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we get the final formula for the total flux across the cell 
boundary: 

F,, = -vri?i/2A(i5Si?i/2). (9) 

Note now that if S and D do not depend on radius than 
Fab = —T^DYi < i.e. the disk becomes an accretion disk. 



2.3. The equation of the difFusion disk evolution 

The law of the mass conservation means: 

dFab 



27ri?AE 



Using Eq. ([9]) we find finally: 



OR 



■At, 



'dt 



R'dR 



(10) 



(11) 



This is the required equation which describes the evo- 
lution of the diffusive disk. 

If we take v — D/Z then we can cast Eq. pT|) into the 
familiar equation for accretion disk evolution (see Pringle 



'dt 



i_d_ 
RdR 



d_ 

dR 



(T,iyVR^ 



(12) 



This equation can be derived from the Navier-Stokes equa- 
tions for viscous disks, where v is the standard viscosity. 
Since D is the true difFusion coefRcient of the matter, we 
have found that the ratio ly/D must be 1/3 in our picture. 

3. The difFusion equation for Keplerian disk: 
multicomponent medium 

Another way to 'derive' the ratio oiv/D is to start from the 
Ansatz that each gas or dust species evolves independently 
according to: 



dt 



RdR 



(13) 



where is the surface density of species i. The E^ obey 
completeness: 



Si 



E„ — E 



(14) 



where n is the total number of species. 

The underlying argument to support the Ansatz of 
Eq. ([Tn)l is the following: Suppose that gas-t-dust parcels 
indeed split up and move inward/outward, as assumed in 
Sect. [21 The relative abundances of the various species are 
the same in the inward moving parcel as in the outward 
moving parcel (they were both from the same original par- 
cel). This means that if the angular momentum stays con- 
stant between the inward and outward moving parcel, the 
same can be said of the individual species. So: species i in 
the inward moving parcel loses as much angular momentum 
as the same species in the outward moving parcel gains. 
Therefore the diffusion happens in each species identically, 
and therefore we argue that Eq. holds for each species 
separately, independent of the other species. This is even 
true if one of the species represents 99.9% of the mass and 
the other only 0.1%. 



So, starting from Eq. p3p we can rewrite 

dY^ _ 3__d_ 
~dr " R'dR 

- 

^ R'dR 



(15) 



If we write Eq. as a conservation equation as E -I- 

{l/R)d[RY.Vfj]/dR — where vji is the average radial ve- 
locity of the matter (including all species) , then by Eq. 
this velocity is given as: 



3 d 



T,iyy/R 



Using Eq. ^B)) in Eq. ([TS]) one obtains: 



(16) 



(17) 



This i s the same as the equation introduced by Morfill & 
Volk (|1984D but now with: 



V 

D 



1/3 



(18) 



So we have shown here that, provided i'/D=l/3, the 
picture proposed by Morfill & Volk (1984 ) that diffusion 
takes place in the comoving frame of an otherwise laminar 
accretion disk is equivalent to our picture of exchanging 
Keplerian parcels. 

Again we stress that our derivation is based on a simpli- 
fying assumption. One can consider our value of v/D = 3 as 
a minimum value of i/fD. At the other extreme is the pic- 
ture of a perfectly laminar but viscous disk with i^/D = inf. 
The true value must be within these limits. 

4. Effect of z^/I? on abundance of crystalline 
silicates in disks 

Although our simple analytic argumentation predicts that 
Sc = 1/3, MHD simulations and experiments with rotat- 
ing turbulent flows show other ratios. At the moment we 
need to consider values ranging from Sc = 1/3 ( this study) 
to Sc = 10 (Carballido, Stone & Pringle l2005t Johansen, 
Klahr & Mee [2006') . So how do these various values affect 
the distribution of molecular species and dust species in the 
disk? What are the consequences of this, for instance, for 
the problem of crystalline silicates in disks? One may argue 
that the difference between v/D—l/S and v/D — 1 is 'only' 
a factor of 3 and therefore not very strong: the time scales 
of diffusion will be different by only a factor of 3. However, 
for accretion disks in semi-steady state, a difference of 3 in 
the diffusion constant can have very strong effects on the 
distribution of thermally processed (crystalline) particles 
in the disk. If we assume that particles are only thermally 
processed in the very inner regions of the disk (i.e. close 
to the star) then to find par ticles at larger radii, the y hav e 
to diffuse outward (e.g. Gaill2001|; Bockelee-Morvan |2002D . 
The efficiency of this outward diffusion is very dependent 
on the v/D ratio. This is due to the fact that the trace par- 
ticles (crystalline dust) have to 'swim upstream' against the 
accretion flow to r each larger radii. This was modeled by 
Clarke & Pringle ()1988[ ) who showed that the abundance 
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of such a tracer as a function of radius can be derived an- 
alytically for steady accretion disks. Since in the absence 
of disk instabilities the inner regions of viscously evolving 
disks change over a time scale much larger than the local 
viscous time, these inner regions (say, inward of a few tens 
of AU at an age of 1 Myr) can be regarded as being in 
semi-steady state. This simplifies the analysis drastically. 
We can then start from the standard advection-diffusion 
equation, after Morfill & Volk 



d /Si 



(19) 



For the special case of v/D=l/2> this becomes equal to 
Eq. (|17p . but we will retain the more general form of 
Eq. (fTOj) here. For a steady state disk we have (for R ^ i?in 
where i?in is the inner radius of the disk): 



M 



with 



^ttRYiVe., 

3 1^ 

' ~2R' 



(20) 



(21) 



Here M is the accretion rate of the disk. With the Schmidt 
number Sc = i^/D, by using Eqs. (|20ll2ip . and by taking 
dUi/dt = (steady state) we can rewrite Eq. ([19]) in the 
form: 

' ' ' ' (22) 



d\gR V S 
The solution is: 

s 



3Sc9(lgi?)2 V E 



.(0) 



.(1) 

R 



(23) 



This is a general steady-state solution for mixing of any 
tracer in the disk for R ^ i?i„. The af^'^ is a background 

abundance while a'"^^ is the abundance at some radius 

In general, the Eq. ([23]) is appropriate only if no large- 
scale radial flows are present in the disk. However, it is 
possible to modify this equation in a way that it approxi- 
mately accounts for the effect of large-scale advective flows 
onto emerg ent spectra. Following t he res ults of Keller & 
Gail (|2004l ) and Tscharnuter & Gail (|2007t ) we assume that 
near the disk's midplane the material moves outward, while 
in the surface layers the material flows inward. Since we as- 
sume that the surface layer of the disk is responsible for 
the emergent spectra we change the Eq. ([^ so as it de- 
scribes the abundances of the tracers in the surface layer. 
Therefore, we substitute Vr{1 + C) for Vr in the Eq. pi)) , 
where C is the circulation velocity in units of the vertically 
averaged accretion velocity, and thus we come to a more 
general relation for the tracer abundances in the surface 
layer: 

X;^ / |sc(i+c) 



Ri 
R 



(24) 



We suppose that C can easily be 2-3 or more which has 
major consequences for the mixing. Given this modification, 
we introduce the effective Schmidt number: 



Sc = Sc(H-C). 



(25) 



Note, that large values of C 'mimic' large Schmidt num- 
bers, Sc, in the picture with no circulation. 




0.0 100.0 1000.0 

R [AU] 

Fig. 3. The abundance of crystalline silicates in a steadily 
accreting protoplanetary disk around a Herbig star of Af, = 
2.5M0, i?* = 2.5L0 and T, = lO-^K. Solid lines: assum- 
ing that only thermal annealing in the warm inner regions 
(here: inward of about 1 AU) can produce crystalline sili- 
cates. Dotted lines: assuming that some other mechanism 
has produced a global base level of crystallinity of 10%. 



Let us now apply these results to the problem of crys- 
tallinity of dust in a disk which is hot enough for thermal 
annealing inward of i?anncai- Since the dependence of the ef- 
ficiency of thermal annealing on temperature is extremely 
sharp, one can say that inward of i?anncai all dust is crys- 
talline. For R > i?anneai there is no annealing and the mix- 
ing solution Eq. (|24p is valid in which i?i = i?annoai and 



CT^Jyst — 1. The background crystallinity <y)^^Lt is the glob 



(0) 



cryst 



cryst 



ally present level of crystallinity, discussed below. One sees 
that the abundance of crystalline silicates at any given ra- 
dius depends only on i?anncai, the base level of crystallinity 
CT^°ygt and the value of Sc. It does not depend on the value 
of V itself (i.e. not on the value of a). For strong or weak 
turbulence the results are the same, as long as Sc remains 
the same. The solutions to Eq. (I24p for a-f^\st — (no global 

crystallinity) and for a^y^t = 0.1 (10% global crystallinity) 
are shown in Fig. [3l 

It is important to keep in mind that these analytic 
solutions for the level of crystallinity are only valid for 
i?i = i?anneai ^ ^in- This is the case for Herbig stars, since 
they are so bright that the annealing radius is around 1 
AU while the inner disk radius is more than hundred times 
smaller. For T Tauri stars, on the other hand, the annealing 
radius is only about ten times larger than the inner radius, 
which means that the effects of the no-friction boundary 
condition of the disk at i?in are still strongly affecting the 
solution, and hence the above analytic solution is not com- 
pletely correct. Nonetheless, the rough form remains correct 
even in that case. 

Now let us focus on the case without a base level of 
crystallinity. It can be seen that the abundance of crystals 
at large radii is extremely sensitive to the power index Sc. 
For Sc = 10, i.e. for large Schmidt numbers or large surface 
inflow velocity, there is essentially no outward mixing. For 
Sc = 1 or Sc = 0.7 the crystals are mixed out to a level of 
3% and 10% respectively at 10 AU. For S"c = 1/3 the 10% 
crystallinity is reached at 100 AU, i.e. out to a hundred 
times larger radius. This shows that the seemingly small 
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differences between the various theoretically derived values 
of Sc make a very large difference in the distribution of 
crystalline silicates in the disk. 

So far we have assumed that crystalline silicates are only 
produced by thermal annealing in the inner disk regions 
(i? < i?anncai)- Howevcr, there is an on going debate about 
what causes dust grains to become crystalline. While t he 
thermal annealing and radial mixing of dust (Gaill2001f) is 
a natural phenomenon that is likely to happen in all disks, 
there may be additional mechanisms o f crys tallizing dust 
such as s hock h eating (Harker & Desch [20021 ) or lightning 
fGibbard ll997l and references therein). Evidence for these 
alternative mechanisms is found in infrared spectroscopy 
of protoplanetary disks, where it is found in many sources 
that the crystalline grains are dominated by forsterite in the 
outer regions of many disks. This is contrary to the radial 
mixing models which predict that enstatite dominates in 
the outer disks (Bouwman et al. in prep.). 

From the solar system there is also some evidence that 
other mechanisms might crystallize dust or dust aggre- 
gates (e.g. review by Trieloff & Palme 2006). However, the 
inner regions of disks, near the dust sublimation radius, 
are always warm enough to assure that thermal anneal- 
ing takes place. Therefore, independent of the existence of 
other crystallization processes, we expect with near cer- 
tainty that inward of some radius all the dust is 100% crys- 
talline (i.e. (Tcryst + '''cryst = !)■ The othcr crystallization 
processes, assuming that they can crystallize dust through- 
out the disk, produce a base level Ccryst of crystallinity. The 
dotted lines in Fig.[3]show how the combined effects of crys- 
tallization work out. We assume here that the time scales 
of these additional crystallizing processes are longer than 
the local viscous time in the region we study here, so that 
the level of crystallinity in the outer disk regions has been 
built up over a long time and remain nearly unchanged on 
the shorter time scales in this inner region. 

Another source of a 'base-level' of crystallinity cTcryst 
that thermal annealing may also happen throughout the 
disk in the v ery ea rly disk-formation phases (DuUemond, 
Apai & Walch [20061 ) . They showed that in the stage that the 
disk is formed, depending on initial conditions, much of the 
outward transport of material may not be due to the 'up- 
stream' diffusion, but due to the fact that the disk materia l 
itself is spreading outward (Lynden-Bell & Pringle [l974l) . 
Later, this pre-processed material accretes back inward and 
forms a base level of crystallinity (Tcryst (see DuUemond et 
al. [20061 . Fig. 2). 

If large-scale circulations are present in the disk then 
over a large distance there may be vertical mixing between 
the outward moving midplane and inward moving surface 
layers. This is likely to produce an enhanced 'base' value of 
crystallinity as it inserts crystalline grains into the surface 
layers at large radii. 

The conclusion of this section is that radial mixing can 
have a strong effect on the distribution of crystalline dust, 
which is important for understanding infrared observations 
of protoplanetary disks, as well as for understanding the 
constituents of meteorites. Radial mixing may also strongly 
affect disk chemistry or grain growth. A small difference in 
the assumed value of Sc as well as of C could make a strong 
difference in the results of the models. 



5. 'Measuring' Sc from disks around Herbig Ae 
stars 

It is clear that it is of high importance to understand which 
is the precise value of Sc. Theories will likely always have 
uncertainties. In this section we propose a way to derive 
the value of Sc ohservationally. From Fig. O it is clear that 
for Herbig Ae stars the crystallinity of the dust in the re- 
gion roughly between 1 AU and 10 AU is strongly affected 

by the ratio Sc, even if there is a base level Ccryst crys- 
tallinity due to other processes. And in the case of Schmidt 
numbers below 1 the radial mixing may even affect regions 
as far as out 100 AU. If we can ohservationally measure the 
crystallinity as a function of radius then we may be able to 
determine both Sc as well as cF^c2/st- 

The abundance of crystalline dust in a disk can be de- 
termined from t he an alysis of infrared spectra (see review 
by Natta et al. |2007[ ) . For ground-based observations the 
8-13 ixm window is well-suited for this. There are by now a 
number of instruments and telescopes that can spectrally 
and spatially resolve disks in the 10 micron regime at spatial 
resolutions of about 20 AU for typical sources: e.g. VISIR 
at the VLT, COMICS at Subaru and T-ReCS at Gemini- 
South. A spatial resolution of about 1 AU can be obtained 
with mid-infrared interferometry with the MIDI instrument 
on the VLT (Leinert et al.'2004). Although in particular the 
interferometry observations do not directly yield spectra as 
a function of radius, they do give information on those spa- 
tial scales that can be compared to model predictions. Of 
course these model predictions have to be done for each 
source individually. 

To demonstrate the principle, we show the intensity- 
spectrum at various radii in a model of a disk around a 
Herbig Ae/Be star of = 2.5M0, i?* = 2.5^© and 

— 10'*K, i.e. the same system as was shown in Fig. [3| 
We took a disk with a mass of 0.0015 M©, a size of 90 AU 
and a surface density profile of S(i?) oc . The 2-D ax- 
isymmetric density structure of the model was taken to be 
p{r,z) = S(i?)exp(-^72i?2)/(^rp(i^)^/2^), where Hp{R) 
is the pressure scale height estimated from the midplane 
temperature of a si mple fl aring disk model of the kind of 
Chiang & Goldreich (|l997f ). We included two dust species, 
both consisting of 0.1 /xm size grains. The first species is a 
mixture of 25% amorphous olivine, 25% amorphous pyrox- 
ene and 50% carbon. The second species consists of 25% 
crystalline enstatite and 25% crystalline forsterite and 50% 
carbon. Within each species the various components are 
in thermal contact, but the two species are not in ther- 
mal contact with each other . We used optical constants 
from Dorschner et al. (|l995l ) for the amorphous olivine 
and pyroxene, from Servoin & Piriou_p,973) for crystalline 
forsterite and from Jager et al. (|1998l ) for crystalline en- 
statite. The method of computation of the opacit ies from 
these constants is described in Min et al. (|2005[ ). Given 
the density structure, the opacities and the stellar param- 
eters we can apply a Monte-Carlo radi ative transfer code 
(RADMC, see DuUemond & Dominik I2004D to compute 
the dust temperature everywhere in the disk. RADMC is 
a Monte Carlo code for dust continuum radia tive tr ansfer 
based on the method of Bjork man & : Wood (|2001l ) com- 
bined with the method of Lucy (|1999[ ) . This gives the typi- 
cal structure of irradiated disks which have a warm surface 
and a cooler interior. A ray-tracing program is then applied 
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Fig. 4. Intensity spectra of the disk (arbitrarily normalized) at different radii, for different assumed Sc, for the standard 

model shown in Fig. [3] Here the base level of crystallinity cr^|!ygj is chosen to be 0, i.e. the spectra correspond to the 
crystallinity levels shown as solid lines in Fig. [31 




to compute, from these density and temperature structures, 
a spectrum. Due to the warm surface layer lying on top of a 
cooler interior, the dust features appear in emission in the 
spectrum. Figure |4] shows the intensity spectra {Iu{R)) ob- 
tained in this way at various radii (arbitrarily normalized) 
for four different values of Sc, assuming that there is no base 
level of crystallinity, i.e. that thermal annealing in the inner 
disk regions is the only source of crystallinity. The intensity 
spectrum is the intensity of the face-on disk image as seen 
at that particular radius, as a function of wavelength. In 

Fig. [5] the same is shown, but including a base level cr^cryst 
of crystallinity due to either the disk's formation history 
(e.g. DuUemond et al. I200 6D or due to additional crystal- 
lization processes (e.g. Marker & Desch 2002). 

From Figs, m and [5] one can see that both the base level 

tr^J^yg^. and the Sc affect the spectra in a way that would 
allow the reconstruction of these two values from the obser- 
vations. The CTcryst is reflected in the fact that the strength 
of the crystalline features decreases with radius but gets 
'stuck' at some level which is set by the value of o^^rst- The 



Sc value is reflected in the rapidity by which the crystalline 
features become weaker with radius as one goes to larger 
radii. And i?anncai is found as the radius by which the crys- 
tallinity suddenly starts to drop. In principle the -Ranncai 
should follow directly from theoretical considerations, as 
the dust temperature in the disk's surface layer, for a given 
grain size, can be determined theoretically. However, the 
grain size is not known beforehand, and an uncertainty of 
the grain size in the 0.1 to 1 /im regime can not be con- 
strained by the analysis of the 10 /im feature shape, but it 
does influence the dust temperature. Therefore an observa- 
tional determination of i?anneai is necessary and possible. 

However, there are some difficulties with this scheme. 
One obvious problem is that one requires rather accurate 
measurements of the spectrum at many radii. Since most 
of the action takes place very close to 1 AU, these mea- 
surements would rely strongly on the MIDI interferometer 
at the VLT or its planned successor, and it would require 
measurements with at least 4 baselines that are all ideally 
aligned in the same position angle on the sky. Moreover it 
would be best if the disk is nearly face-on and bright. 
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Another complication is that thermal annealing may 
take place in the surface layers at a different radius than 
the midplane. If the disk's heating is dominated by irradi- 
ation, then the midplane temperature of the disk is much 
smaller than the surface temperature. In this case there 
are radii where the dust is above annealing temperature in 
the surface, but below annealing temperature in the disk's 
midplane. It then depends strongly on the vertical mix- 
ing efficiency whether the annealing in the surface can also 
strongly affect the crystallinity of dust in the midplane. If 
not, then the i?anncai relevant for the radial mixing might be 
smaller than the radius where the surface is annealed. One 
would then expect, as a function of radius, the crystallinity 
to be unity out to that surface annealing radius and then 
to suddenly drop to the level of the power-law which origi- 
nated at a smaller annealing radius. In principle this could 
also be observed, if sufficiently finely-spaced measurements 
are taken. If accretional (viscous) heating could heat the 
disk's midplane sufficiently strongly, then the oppose may 
be true: the annealing radius would then be larger than 
the annealing radius computed with irradiation only, and 
only drop at larger radii. Moreover, if meridional flows are 
present, the picture might blur even more. 

It is clear that the full interpretation of these measure- 
ments will not be easy. But since the Schmidt number is 
such a fundamental number for accretion disk theory, it 
might be worthwhile to attempt to do such measurements 
nonetheless. 



6. Conclusions 

In this paper we discuss the ratio between turbulent dif- 
fusion and viscosity in protoplanetary disks and study an 
effect of this ratio on the dust crystallinity. Our conclusions 
can be summarized as follows: 

1. The efficiency of radial mixing of any tracers in a 
steady state disk depends on the ratio between viscos- 
ity and diffusion coefficients i.e. on the Schmidt number 
Sc = i^/Z) if we assume that the disk is well mixed in 
vertical direction. Therefore, the Schmidt number plays 
an essential role in the understanding of the distribution 
of the various chemical and dust species in disks. 

2. We present the simple analytic model of the diffusively 
spreading disk. We show that evolution of the diffusively 
spreading disk and evolution of the viscousely spreading 
disk are described by the same equations if we assume 
iz/D — 1/3 where f is the effective viscos ity coefficient 
in the classical Shakura & Sunyaev lUll) picture. This 
does not mean, however, that we necessarily argue that 
true accretion disks have this Sc = 1/3: this depends 
on the validity of our assumptions (vertically averaged 
disk, good mixing between dust and g particular 
kind of turbulence). However, one can regard the value 
Sc = 1/3 as the lower limit of the Schmidt number 
which can be included in the models of radial mixing in 
Keplerian disks. 

3. If large-scale meridional flows are present in the disks 
then the Schmidt number cannot be used as a direct 
measure of the mixing efficiency. However, it is possible 
to approximately describe the abundance of the trac- 
ers in the surface layers of the disk by introducing the 
effective Schmidt number Sc = Sc(l -I- C) where C is 



the layer velocity in units of vertically averaged accre- 
tion velocity. This approximation is valid as long as no 
vertical mixing between in- and outward moving layers 
is assumed. In this case, we expect less abundance of 
crystalline dust in the surfacejayers. 
4. We show that the value of Sc may possibly be mea- 
sured in protoplanetary disks around young stars by 
using mid-infrared interferometry. This requires mea- 
surements of a single source at many baselines which all 
lie along the same position angle across the sky. 
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Appendix Al 

Consider Eq. ^ , and produce the Taylor expansion around 
central cell, R: 

MRV = Ma{R + h)(v+^h+ \^h^ 



(26) 



After elementary operations, using Eq. ([1} and dropping 
terms of order ft,"^ we get: 

,„,-M„).^ii|i^^ (27) 

Si""'' 



For the Keplerian disk, V ^ R this relation becomes: 

(28) 

Using Eq. ([1]) we obtain: 

M..f(:-ii). 



Appendix A2 

Rewrite the fluxes M^^^ and A/f^+''^ using Eq. © and (gl): 



M(«) = iM(«)(l + i- 
2 \ 4i? 



(31) 



(32) 



where M(^) and A/(^+'') are the masses emerging from the 
cells with radii R and R + h respectively. Denote Fab as a 
total flux across the boundary between the cells: 



(33) 



Ya. Pavlyuchenkov and C. P. DuUemond: Dust crystallinity in protoplanetary disks 



9 



If we take the Taylor expansion of Eq. around R and 
leave only the linear terms in h we get: 

^ ^ h flM dM\ 



It is convenient to rewrite the last equation in the form: 

(35) 



Appendix A3 

Let us comment on an important point in the above deriva- 
tion. One may argue that the factor 1/2 in Eq. ([5]) is an 
arbitrary number. Thus, the factor which appears in the 
right part of final Eq. (fTTjl (now it is 1) is also arbitrary. 
However, we support our choice of 1/2 in Eq. ([H]) based on 
the following arguments. It is easy to rewrite the derivation 
for ID planar geometry i.e. to obtain the diffusion equation 
in Cartesian coordinates. In terms of the above equations 
we should simply set Ma — Mb = M/2 instead of using 
Eqs. ©-(g]). The Eq. jS]) is replaced by 



hdM 
2 aR ' 



and Eq. ([7]) is rewritten as 

M = S • K 



turb 



•At, 



(36) 



(37) 



where M is the total mass of the current cell. We should 
substitute these relations into analogue of Eq. pU]) for pla- 
nar geometry 

Ai? • AS = -{FabiR + h)- FabiR))At. (38) 

In this way we obtain the classical equation: 

^-D^ (39) 
dt~^dR^' ^^^^ 

only if we introduce the diffusion coefficient, D, in the form 
of Eq. ([5]). In other words, we choose the convention Eq. ([5]) 
in such a way that Eq. pT|) turns into Eq. p9|) in the limit 
of i? ^ oo and _D=const. 
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